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Variation of magnetocaloric (MC) effects has been systematically investigated for colossal mag- 
netoresistive (CMR) manganites _Ro.6Sro.4Mn03 (7?=La-Gd). As the one-electron bandwidth is 
reduced, the temperature profile of MC effect, i.e. field-induced entropy change, exhibits a steeper 
drop below the ferromagnetic transition temperature due to its first-order nature promoted by a 
competing charge-orbital ordering instability. For these small-bandwidth systems adjacent to the 
metal-insulator phase boundary, a rectangular-shaped profile for the entropy change emerges with 
an anomalously wide temperature range and a considerable magnitude. Model calculations have 
indicated that the fluctuation enhanced in the phase-competing region has a strong impact on such 
MC features, which can be extensively controlled by the chemical composition. 



A magnetocaloric (MC) effect refers to the isother- 
mal entropy change, AS", and the adiabatic tempera- 
ture change, AT a( j, induced by applying (or removing) 
a magnetic field to the materials. Since it constitutes 
the basis of the magnetic refrigeration with environmen- 
tal friendliness and high efficiency, the quest for giant 
MC materials is currently of great interest i It was re- 
cently found that systems undergoing a first-order mag- 
netic transition show a large AS at the transition tem- 
perature, Tq . Alloyed metals with Tq near room temper- 
ature were hence the main subject of research, aiming at 
the application to refrigerators and air-conditionersi^r— 
In addition, the low-temperature refrigeration such as 
hydrogen liquefaction now requires the wide-range refrig- 
eration ability down to ^20 K. In such contexts, colossal 
magnetoresistive (CMR) manganites, where Tq of ferro- 
magnetic metallic (FM) phase can be widely tuned by 
changing the one-electron bandwidth, can be promising 
MC materials. Although various large-bandwidth man- 
ganites have been intensively investigated so far^ system- 
atic studies focusing on the narrow-bandwidth systems 
with low Tq are quite rare££ In the latter case, as a re- 
sult of phase competition with a charge-orbital-ordered 
(CO/ 00) insulator, sharp FM transitions are observed 
as the CMR phenomena^ where substantial AS" may be 
generated by a weak field. 

The purpose of this Letter is to reveal the compre- 
hensive feature and the underlying new physics for the 
MC effect of the CMR manganites. As an ideal arena 
for this, we have adopted a series of single crystals 
of i?o.6Sro.4Mn03 (i?=La-Gd), which exhibit the wide- 
ranging Tq from 375 K to 59 K by decreasing the ionic 
radius of R site, or the bandwidth [Fig. |3fa)]. The mag- 
nitude of the bandwidth is determined by a tolerance 
factor /, defined such that /= (rA + ro)/\/2(i'Mii + ''o). 
with 7"a, fMni and ro being the (averaged) ionic radii 
of the perovskite A- and Mn-site cations and oxygen, 
respectively.^^ Even in the most reduced-bandwidth 
regime, no long-range CO/OO phase is realized due 
to the quenched-disorder effect arising from large mis- 



match in ionic size of R and Sr, and instead a spin 
glass (SG) phase shows up with only short-range CO/OO 
correlation^ We have revealed a systematic change in 
MC properties throughout the phase diagram and, in 
particular, an anomalous MC effect characterized by a 
wide temperature range exhibiting substantial AS*. De- 
tailed comparison with model calculations indicates that 
the origin is ascribed to the gigantic phase fluctuation 
evolving near the phase boundary to the SG phase. 

Single crystals of i?o.6Sro.4Mn03 (i?=La-Gd) were 
grown by a floating zone methodi 12 i 13 Powder x-ray 
diffraction exhibited that the obtained crystals are of 
single phase. The magnetization was measured with a 
superconducting quantum interference device (Quantum 
Design). The specific heat measurement was performed 
by a standard heat relaxation method, using Physical 
Property Measurement System (Quantum Design). 

We first focus on the i?=Sm compound, which ex- 
hibits a first-order FM transition around 123 K [Fig. 
[TJc)], as a typical example. We can indirectly ob- 
tain AS by utilizing the Maxwell relation: AS(H, T) = 

f"(dM/dT) H dH, where H, T, and M are applied 
field, temperature, and magnetization, respectively. Us- 
ing all the discrete M-H data shown in Fig. QJa), we 
numerically calculated the AS values [open circles in 
Fig. [He)]. As an alternative, measurements of the spe- 
cific heat C(H,T) at H=0 T and 6 T [Fig. QJb)] can 
also figure out AS based on the thermodynamic law: 
AS(H,T) = J^[C(H,T)-C(0,T)}/TdT. Because of the 
difficulty in measuring the sharp peak stemming from the 
latent-heat contribution at T by the relaxation method, 
we interpolated the data by a Lorentzian function fol- 
lowing the literature^ when performing the numerical 
T integration. Thus calculated AS values [closed cir- 
cles in Fig. [TJc)] are in good agreement with those from 
the M data. This consistency between the magnetic and 
calorimetric measurements ensures the validity of relying 
on the Maxwell relation in the present systems, despite 
their first-order nature of the transition^ In the follow- 
ing, we have estimated AS from the M data. 
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FIG. 1: (Color online) Magnetocaloric properties for R—Sm. 
(a) Isothermal magnetization (M) versus magnetic field (H) 
up to around To, taken with the temperature (T) interval of 
2 K (not all shown here) in the H- and T-decreasing runs, (b) 
Temperature profile of specific heat (C) at T and 6 T. The 
dashed line is the Lorentzian fitting to the latent-heat peak, 
(c) Temperature dependence of entropy change for H—6 T, 
calculated from both M and C data. The corresponding M 
at 0.5 T is also shown. 



We show the temperature dependence of M for 
i?o.6Sro. 4 Mn0 3 (i?=La-Gd) at 0.5 T in Fig. Ufa). For 
the i?=La crystal with the largest bandwidth, a sim- 
ple second-order FM transition was discerned at Tq as 
high as 375 K. With decreasing the bandwidth, Tc sys- 
tematically decreases and the FM transition becomes 
steeper. Judging from the temperature hysteresis at the 
transition, the critical end point of the first-order transi 
tion seems to locate in the vicinity of i?=Smo.sNdo.5 
In the systems with smaller /, such as R=G&i- x Syh x 
(0.5 < x < 1), the apparent first-order FM transition was 
observed, except for the R=Gd crystal (SG insulator). 

The resultant curves of —AS" are shown in Figs. &b) 
and (c) for H=l T and 6 T, respectively. The tempera- 
ture range of the significant MC effect can be widely con- 
trolled from 50 K to 400 K by changing the R species. 
In the large-bandwidth systems (i?=La and Pr), — AS 
exhibits a broad peak at Tq. Both its height and width 
steadily increase with increasing the magnitude of H, 
characteristic of the second-order transition. When the 
bandwidth of the system is reduced, the shape of — AS 
becomes sharper with a sudden drop below Tq. Accord- 
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FIG. 2: (Color online) Temperature profiles of M at 0.5 T 
(a) and AS for H=l T (b) and 6 T (c) for the single crystals 
of J?o.6Sro.4Mn03 (i?=La-Gd). The solid and dashed lines 
in (a) correspond to warming and cooling runs, respectively. 
The AS data was obtained in the H and T-decreasing runs. 



ingly, the peak value becomes larger and reaches ~10 
J/kg K (~18% of the total magnetic entropy) for H=6 
T for i?=Sm, although it saturates for i?=Gdi_ a; Sm a ; 
(0.5 < x < 1). There, the height of the -AS peak is 
not sensitive to H above 1 T, whereas its width increases 
nearly linearly with H, as is typical of the first-order 
transition. Noteworthy is that their peak shape changes 
to the rectangular one with an anomalously wide plateau 
temperature region (e.g., ^50 K for i?=Gdo..5Smo.5 for 
H=6 T), as compared with the case of rare-earth alloys 
showing the first-order transition at similar Tq ^ A grad- 
ual but large increase in — AS* below 40 K is probably due 
to the contribution from the Gd spins. 

To provide a perspective of the MC properties for 
i?o.6Sro.4Mn03 (i£=La-Gd), we display in Fig. [3jb) a 
contour plot of —AS" for H=6 T on the entire f-T plane. 
Contour values arc determined by interpolation based on 
the data shown in Fig. HJc). Depending on the behavior 
of AS, we roughly classify their MC effect into three re- 
gions in the f-T plane. In the region (I) corresponding to 
i?=La-Pr with the second-order FM transition, —AS" ex- 
hibits a broad peak around Tq, forming a yellow-colored 
(online) hump. In the intermediate-bandwidth region 
(II), such as i?=Nd-Sm .5Ndo.5, Tq begins to rapidly 
decrease with decreasing the bandwidth due to the en- 
hanced antifcrromagnctic CO/OO correlation^ Further- 
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FIG. 3: (Color online) (a) Phase diagram for 7?o.6Sro.4Mn03 
(i?=La-Gd) as a function of a tolerance factor /, which is 
defined in the text. The ferromagnetic metal and spin-glass 
insulator are denoted by FM and SG, respectively. The ferro- 
magnetic transition temperature 7c, determined as the peak 
position of —AS, is represented by closed circles, while the 
SG transition temperature T g by closed triangles, (b) Con- 
tour plot of -AS for H=6 T on the f-T plane. 



more, the FM transition approaches the first-order one 
accompanied by a sharp drop in —AS* below Tq, which 
results in the evolution of the red-colored (online) sharp 
ridge. In the region (III) corresponding to i?=Gdi_ a; Sm a ; 
(0.5 < x < 1) with the narrowest-bandwidth, Tc is 
markedly suppressed as the SG phase is approached. 
These crystals undergo a first-order transition even at 
high fields, leading to a steep rise in —AS above T> 
The resultant rectangular-shaped —AS* is highlighted as 
the red-colored (online) plateau structure below ~1 50 K. 
In each case, we hereinafter make a detailed comparison 
between the model calculation and the observed data. 

In the region (I), the mean-field calculation of AS using 
the Heisenberg model of spin s=l .8 with the equivalent 
2c can almost reproduce the experimental values, as ex- 
emplified for i?=Pr in Fig. 0|a). Because of the large 
spin value and the weak CO/OO instability, the mean- 
field approximation works well in this regime. Neverthe- 
less, a small deviation is discerned particularly above Tc, 
which can be improved by the theory taking account of 
the bicritical fluctuation (vide infra). 

As shown for iJ=Smo.sNdo.5 in Fig. Hfb), the Heisen- 
berg model can no longer reproduce the sharp drop in 
—AS below Tc in the region (II). To compare this be- 



havior with the first-order transition model, we took 
account of the magnetostriction effect in the mean- 
field Heisenberg HamiltonianJ^— Hmf = [~2zJ(£)(s) — 
gu B H]j: iSl + j: {h]) [2J(e)( S ) 2 + (n/2)(l-e) 2 }, where e 
is the amount of the dimcnsionlcss lattice strain, z the 
number of the nearest neighbors, (s) the thermodynamic 
average of the spin, Sj the spin operator at the i site, 
g the Lande factor, ub the Bohr magneton, and k the 
elastic constant. The magnetostrictive exchange energy 
is given by J(e) = Jo — Ae, where Jo is determined as 
2zJq =3A;b@p/s(s + 1) (fee is the Boltzmann constant) 
and A is a positive coefficient. Op corresponds to the FM 
transition temperature, provided that the one-electron 
bandwidth is simply reduced without the CO/OO insta- 
bility [see Fig. EJa)]. The elastic energy term f (1 — e) 2 
effectively includes the sum of the energy loss due to the 
lattice strain and the energy gain due to the Jahn- Teller 
effect, arising from the CO/OO correlation. Because 
e = 1 (for H=0 T) above Tc so as to minimize this term, 
the reduced ferromagnetic exchange coupling Jo — Ae re- 
sults in Tc lower than Op , which mimics the short-range 
CO/OO correlation evolving above Tc- After calculating 
(s) and e by minimizing the free energy, we determined 
A and K so that Tc and the temperature profile of AS 
are most consistent with the experimental results. The 
result fitted for i?=Smo.sNdo.5 is shown by a dashed-two- 
dotted line in Fig. SJb). Whereas it reproduces the —AS 
drop below Tc, its peak height is ^2 times as high as the 
experimental value. Furthermore, the predicted behav- 
ior above Tc differs significantly from the experimental 
results. 

To correctly take into consideration the effect of fluc- 
tuations enhanced in the phase-competing region (II) 
beyond the mean- field scheme, we now deal with the 
phcnomcnological Ginzburg-Landau model based on the 
renormalization group analysis.— This method not only 
successfully explains the first-order FM transition around 
the bicritical point for CMR manganitcs, but also de- 
rives the universal scaling relation between H/M and M 2 
(Arrot plot), expressed as H/M = tm + 2g P AiN p f (r p + 
g pM M 2 ), where r M , r p , N p , and g pM are scaling pa- 
rameters, p the order parameter of charge ordering, and 
f'(x) = f(lna; + l). Following the literature^ we per- 
formed the scaling plot of the M-H data for i?=Pr- 
Sm .5Ndo.5 by regarding / as a parameter controlling 
the bandwidth. Note here that tm and r p are ex- 
panded with respect to T and / near the bicritical point: 
fi = CirAT + CifAf (i = M, p), where and Cif are 
constants independent of R. The results of the scaling 
plot for i?=Pr-Smo.5Ndo.5 are shown in the inset to Fig. 
|4jb). All the M-H data fall into one curve with good 
accuracy, reflecting that the bicritical fluctuation is en- 
hanced in this regime. Using the obtained parameters, 
we calculated AS for J?=Smo. 5 Ndo.5 [Fig. gfb)]. The 
result well reproduces the broad tail above Tc as well as 
the steep drop at Tc, which the mean- field calculation 
fails. The enhanced fluctuation makes the system sensi- 
tive to a magnetic field, which may effectively broaden 
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FIG. 4: (Color online) Temperature profiles of AS around Tc 
for H=6 T for i?=Pr (a), Sm . 5 Ndo.5 (b), and Gd . 5 Smo.5 (c). 
Results of the mean-field (MF) calculation using the Heisen- 
berg (HB) model and that with magnetostriction effects (MS) 
are shown by dashed-dotted and dashed-two-dotted lines, re- 
spectively, while those calculated using the Ginzburg-Landau 
(GL) theory based on the renormalization group (RG) anal- 
ysis shown by solid lines. Experimental data correspond to 
open circles. Inset to (b) exhibits the scaling plot for i?=Pr 
and Smo.sNdo.s with the curve of f'(x) = |(lna; + l). 

the peak width of —AS. Also for i?=Pr, the renormaliza- 
tion group analysis works better than the mean-field one 
[Fig. Ufa)], indicating some influence of the fluctuation 
even in the region (I). 

The plateau structure of — AS observed in the region 
(III) cannot be reproduced by the mean-field calculation 
in either the Heisenberg model nor that with magne- 
tostriction, as shown for i?=Gdo.5Sm .5 in Fig. life). (We 



are not concerned here with the broad peak at the lowest 
temperature, stemming from the Gd spins.) Although 
the latter model can reproduce a rectangular shape, its 
width is less than 1/3 of that observed and its height 
is about 3 times. This anomalously wide plateau might 
be also caused by the fluctuation further developed in 
the low- 7c systems. As described in the literature^ 
however, the scaling relation cannot be applied to the 
strong-fluctuation regime accompanied by the distinct 
first-order transition like the region (III) . In this regime, 
quenched disorder arising from the A-site solid solution 
completely suppresses the long-range CO/OO phase. It 
leads to the short-range CO/OO state evolving above 
Tcr^ which tends to promote the first-order nature of 
the FM transition. To correctly deal with their MC phe- 
nomena, we hence need a more elaborate model which 
takes account of the strong quenched disorder in addi- 
tion to the phase competition. 

From a practical viewpoint, this series of manganites 
may be suitable as a magnetic refrigerant especially be- 
low ^200 K. There, they indicate the good performance: 
a refrigeration range of as wide as ~50 K with substan- 
tial heat extraction (-AS ~ 10 J/kg K) for H=6 T, 
achieving the cooling capacity-^ of ~500 J/kg compa- 
rable with that of the giant MC materials. Furthermore, 
the plateau-like shape of the —AS peak is suitable for 
the Ericsson cycle, where the constant AS as a function 
of temperature is required within the cooling rcgioni 23 i 24 

In conclusion, we have investigated the overall fea- 
ture of the MC effects for a series of CMR manganites 
-Ro.6Sro.4Mn03 (i?=La-Gd) with the controlled band- 
width. With decreasing the bandwidth (or Tc), the 
— AS peak becomes larger, changing to a rectangular 
shape, due to the first-order nature of the transition in- 
duced by the competition with the CO/OO instability. 
By performing various model calculations, we have at- 
tributed the observed unusually wide temperature range 
of their MC effect to the gigantic fluctuation in the phase- 
competing regime. Thus, the electron correlation can 
strongly affect (or sometimes improve) the MC features, 
giving a new clue to exploring giant MC materials. 

We thank Y. Onosc and H. Katsura for fruitful discus- 
sions. This study was partly supported by Grant-in- Aid 
from Special Postdoctoral Researchers Program and In- 
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